clear all;
t = 0:0.01:16;
e0 = [2 0]';
[t,e] = ode45(@fun0,t,e0);
figure(1);
plot(e(:,1),e(:,2));grid;
xlabel('$e$','interpreter','latex','FontSize', 18);     %注意latex在matlab里的使用
ylabel('$\dot e$','interpreter','latex','FontSize', 18);
figure(2);
subplot(2,1,1);plot(t,2-e(:,1));grid;
ylabel('$c(t)$','interpreter','latex','FontSize', 18);
xlabel('$t$','interpreter','latex','FontSize', 18);
subplot(2,1,2);plot(t,e(:,1));grid;
ylabel('$e(t)$','interpreter','latex','FontSize', 18);
xlabel('$t$','interpreter','latex','FontSize', 18);

%定义函数（非线性分段转化为线性方程）
function de=fun0(t,e)
de1 = e(2);
if(e(1)<-0.2)
    de2 = 0.2*4-e(2);
elseif (abs(e(1))<0.2)
    de2 = -4*e(1)-e(2);
else de2 = -0.2*4-e(2);
end
    de = [de1 ;de2];
end